Structure function of passive scalars in two-dimensional turbulence 
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The structure function of a scalar #(x, t), passively advected in a two-dimensional turbulent flow 
u(x, t), is discussed by means of the fractal dimension 5^ of the passive scalar graph. A relation 



between the scaling exponent of the scalar structure function Z)j" ; (r), and the structure 
function D-2(r) of the underlying flow field is derived. Different from the 3-d case, the 2-d structure 
function also depends on an additional parameter, characteristic of the driving of the passive scalar. 
In the enstrophy inertial subrange a mean field approximation for the velocity structure function 
gives a scaling of the passive scalar graph with Sg < 2 for intermediate and large values of the 
Prandtl number Pr. In the energy inertial subrange a model for the energy spectrum and thus D2(r) 
gives a passive scalar graph scaling with exponent Sg 1 ^ = |. Finally, we discuss an application to 
recent observations of scalar dispersion in non-universal 2-d flows. 
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I. INTRODUCTION 

The dynamics of a scalar field 6*(x, t) advected in a 
turbulent velocity field u(x, t) is of practical relevance in 
many fields of current research such as air pollution or 
chemical reactions in the stratosphere in connection with 
the ozone hole JlJ . Especially for problems in atmospheric 
physics, models of two-dimensional turbulent flows give 
a good approximation of the dynamical processes and 
are frequently used . More recently, two-dimensional 
turbulence has become experimentally accessible in mer- 
cury layers Q , thin salt water layers [||-[| , and soap films 
P-|r2| . Two-dimensional turbulence is also interesting be- 
cause of its fundamentally different behavior compared to 
the three-dimensional case. Since the enstrophy is a sec- 
ond inviscid invariant beside the energy two cascades de- 
velop: starting from a fixed, intermediate injection scale, 
energy is transported to larger spatial scales in an in- 
verse en ergy cascade and to smaller ones in an enstrophy 
cascade @Q- 

The scaling behavior of a passive scalar in a turbu- 
lent fluid was analyzed mainly in three dimensions where 
three different regimes could be identified. Depending on 
the Reynolds number of the underlying fluid turbulence 
and the ratio of the kinematic viscosity to the scalar dif- 
fusivity one distinguishes the viscous-convective Batche- 
lor regime , the inertial-convective regime Jl(],[l7j , and 
the inertial-diffusive regime. In 2-d the situation is more 
complicated, since already the velocity field shows a vari- 
ety of scaling regimes. In particular, the inverse cascade 
process gives rise to the formation of large scale vortices 
that change on very slow time scales only fl8|| and can 
dominate the dynamics of the passive scalar, at least on 
intermediate time scales (l|]2(]] . The formation of coher- 
ent vortices can be suppressed by a large scale dissipa- 
tion mechanism. If this additional dissipation is present 
a statistically stationary homogeneous and isotropic tur- 
bulent flow field develops, that can be characterized by 
its structure function. We assume that a passive scalar 
in such a flow field also develops a statistically station- 
ary state which can be characterized by its own structure 



function. 

The approach used to analyse the structure function 
of the passive scalar is geometric measure theory [^l] . 
This powerful method allows to connect the structure 
function of the passive scalar to that of the underly- 
ing flow field and thus to link the statistical behavior of 
both. The result are scale resolved bounds on the scaling 
behavior. Upper bounds are easiest to derive and often 
give very good results, see e.g. the favorable compari- 
son between theory and numerical simulations in [p5| . 
The derivation of lower bounds is possible |^3| but much 
more difficult and will not be attempted here. So assum- 
ing the reliability of the upper bounds we would like to 
see how the different regimes in u are reflected in the 
scaling properties of the scalar field passively advected 
by the flow. Some aspects of the 2-d case have been dis- 
cussed previously |24], see below. In addition, we would 
like to compare the predictions to the results of experi- 
ments of Cardoso et al. || , where certain discrepancies to 
theory were noted. As we will see the discrepancies can 
be accounted for if the experimentally measured struc- 
ture function is substituted for the velocity field. 

The model we consider is that of a scalar field #(x, i) 
transported in the turbulent flow field u(x, t) according 
to 



r)9 

at 



(1) 



k denotes the diffusivity. The force density fg models ex- 
ternal boundary conditions and the driving and assures 
a statistically stationary field 8(x.,t). The scalar 8 is as- 
sumed to be passive, i.e. it does not affect the dynamics 
and the statistical properties of the velocity field. We 
assume that in the presence of a large scale dissipation 
mechanism a homogeneous, isotropic, and stationary tur- 
bulent state develops. The ratio of the kinematic viscos- 
ity v to the scalar diffusivity n defines the Prandtl num- 
ber Pr = vjn (this is the nomenclature used when 8 is 
a temperature field; if it describes a concentration then 
the corresponding ratio is known as the Schmidt num- 

lg\ 

ber) . The scaling exponents Q of the n-th order scalar 
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structure functions, defined as 



Dl l e \r) = (\6( X + r 7 t)-6( X7 tT)~A> , (2) 

can be obtained from an analysis of the fractal dimen- 
sion 5^ of d-dimensional scalar field graphs; (•) denotes 
the statistical ensemble average. The fundamentals of 
the geometric measure theory approach were laid out 
by Constantin et al. pf]-^] who derived the fractal di- 
mension S g (d is the space dimension). Closely related 
to the present investigation is the application to two- 
dimensional chaotic surface waves [p4| . The Pr depen- 
dence of 3-d passive scalar advection within this approach 
was discussed in [p5[ . As in that work we will aim at a 
rather direct relation between scaling exponents and ve- 
locity structure functions. 

The outline of the paper is as follows. In Sec. |l| the 
basic concepts of the evaluation of the fractal graph di- 
mension are summarized. The results of the mean-field 
approach |2(f | for fully developed two-dimensional turbu- 
lence in the direct enstrophy cascade range - the scal- 
ing behavior of the 2nd order velocity structure function 
D%{r) = (|u(x+r, t) — u(x, t)\ 2 ) - are recalled. In a second 
step we interpolate the scaling of D2 (r) to the inverse en- 
ergy cascade range, where no analytical result is known. 
We obtain -D2M from Fourier transform of an energy 
spec trum as is found in many numerical simulations. In 
Sec. Ill the fractal dimension of the passive scalar graph 



is derived over a broad range of Prandtl numbers, both 
in the enstrophy inertial subrange (ISR) and in the en- 
ergy ISR with the previous relations for the structure 
function. We conclude with a summary, a discussion of 
the relation to the findings in the quasi-two-dimensional 
dispersion experiments by Cardoso et al. ||, and some 
remarks on open questions. 



II. BASIC CONCEPTS 

A. Fractal dimension of the passive scalar graph 

From now on all considerations are made for the case 
of a two-dimensional flow field. The graph of the scalar 



field is then a 2-d surface in 3-d space. The Hausdorff 
dimension of this graph is obtained from the scaling be- 

(2) 

havior of the Hausdorff volume H(G(Br )) of the graph 



G(B^) = {(x,( 
radius r (the 2-d ball Br 

H(G{B™)) 



|x € B {2 \ 6 = 0(x)} over a disk of 
?( 2 h fv5tl 



S (2) 
f 9 



(3) 



(2) 

In two dimensions the fractal dimension S a is connected 



to the scaling exponent C-f\ cf. eq. 
equality [El 



through the in- 



4 2) < 3 - 



C[ e) ■ (4) 

We assume equality in M) |2^,^] and use the relation 
5g = 5 g 2 ^ — 1, where 6g is the fractal dimension of 
the level sets 80 = 6*(x). The relative Hausdorff volume 
H{G(Bi 2] 
ory |2^,|29[ as 

h(g{b1 2) )) 



(2") (2) 

H{GiBr'))/V{B K /') is given by geometric measure the 
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V{B (2) ) 



V{B (2) ) Jb$ 



l + r 2 |V0| 2 d 2 x 
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|V0| 2 d 2 x 



(5) 



(2) 

where the Cauchy-Schwartz inequality and V(B r ) — 
irr 2 were used in the last line. The passive scala r field 
0{x,t) is measured in units of 8 r . m .s. = \J (^ 2 )j thus 
leading to dimensionless 9 = 6/6 r , m , s .. Equation (|5|) is 
a generalization of the well-known volume formula V = 
J ^/gd 2 y to fractal sets, where V is a two-dimensional 
curved hyper surface embedded in the three-dimensional 
Euclidean space and g the determinant of the metric ten- 
sor gtjiy 1 ^ 2 ). 

We now turn to the evaluation of S g 2 \ The term |V6>| 2 
can be replaced by means of (Q) by 



|V0| 2 = — (kA0 2 - (u • V)<9 2 ) 1 fed 
2k 



n9 r 



(6) 



With this eq. (||) becomes 



H(G(B r 2) )) 



< 



V{B {2) ) " \ 




i-[-(u.V)0 2 + K A0 2 ] + -^ 

£j r\j r\i\J ip i> 



d 2 x 



(7) 



We will consider the three integrals under the square 
root separately and denote them by 1%, I2, and ^3, re- 
spectively. In the three-dimensional case [^5| the terms I2 
and ^3 vanish in the large Reynolds number limit. They 
also satisfy the inequality I2 < which changes to 



mated as 



in the two-dimensional case. Is can be esti- 



h = 



,2 r 2 e e 
■a x — 



k9 2 

^r.m.s. 



"jf^Prr 2 , (8) 

r.m.s. 



where the scalar dissipation rate eg = k(|V0| 2 ), the en- 
strophy dissipation rate e u = i/(|Vo;| 2 ), and stationar- 
ity are used. In the case of a three-dimensional passive 
scalar this term contains a factor v 1 ! 2 and thus can be 



2 



neglected. In 2-d the smallest scales are given by the en- 
strophy dissipation rate and this factor disappears. Hence 
I3 cannot be neglected; its importance is evidently con- 
trolled by Prandtl number Pr, length scale r and dimcn- 
sionless prefactor 



-1/3 



92 

r.m.s. 



(9) 



The term I 2 can still be neglected on account of its sub- 
dominant scaling in r. Wc introduce dimcnsionless length 
scales f — r/rj u by means of the enstrophy dissipation 



length 77^ 



,1/2-1/6 



since in 2-d turbulence it is the 



enstrophy cascade that brings the energy to the smallest 
scales where viscosity dominates. 

It follows from (0) for I\ by applying the Gauss Theo- 
rem and the Cauchy-Schwarz inequality 



li < r 

K 



D,,(r) 



(11) 



Combining (||), (0), (||), and we end up with an in- 



equality for the fractal dimension S g 
graph in two dimensions, 



(2) 



of the passive scalar 



2 < 




e 2 

r.m.s. 



Pr f 2 + V~3Pr 



(12) 



where Dm = DJ(el /3 r)l 



This inequality, relating the 

scaling exponent 8)f J to the longitudinal structure func- 
tion of the underlying turbulent flow field Dy is the main 
result of this section. For most of the discussion that 
follows we will assume equality in Jl2]); in the three- 
dimensional case this is a very good assumption |25| . 




dr. (10) 



The quantity u r — 2irr is the circumference. It is possible 
to add Uo, the velocity at the center of Br , due to the 
assumed homogeneity. 

The first term on the right hand side contains the 
square root of the passive scalar flatness. Since we are 
interested in the scaling properties of Ii, it suffices to 
know that the scalar flatness is a constant, independent 
of r. However, there do not seem to be numerical or ex- 
perimental data for the passive scalar flatness in 2-d. 
Data for the velocity field from the experiments |1 and 
the numerical simulations p0[ | suggest Gaussian behav- 
ior in the absence of coherent structures in the regime 
of the inverse cascade. More recent experiments suggest 
that this result also extends into the region of the direct 
enstrophy cascade |3l[]. However, since there are models 
where a Gaussian statistics for a random velocity field 
causes non-Gaussian scalar statistics |p2| , p3[ , this infor- 
mation is insufficient to infer Gaussian statistics for the 
passive scalar. In the following we will work with the 
Gaussian flatness value of three for the passive scalar. It 
should be kept in mind that deviations from this value 
will most likely be scale dependent and will give rise to 
modifications of the scaling exponents. 

The second term is the longitudinal velocity structure 
function Dy(r). Thus we find 



D 2 (f) 



D^(oo) 



Structure functions in two-dimensional 
turbulence 



B. 



To evaluate (|12|) we need information on the scaling 
behavior of the 2nd order longitudinal structure func- 
tion Dy . The longitudinal structure function Dy (r) and 
transversal structure function D± (r) make up the veloc- 
ity structure function D 2 (r) and are connected by incom- 
pressibility, D± = Dm + r-^r. Eliminating the transversal 
part then gives 34 3j| 



D\\(r) 



1 



pD 2 (p)dp 



(13) 



As there are two inertial ranges with several different 
scaling regimes, there is no analytical expression for the 
structure function. As far as we are aware, the best that 
can be achieved analytically is the structure function for 
the enstrophy cascade as discussed by Grossmann and 
Mertens pofl . They used a mean field type approach for 
the fully developed, turbulent velocity field in the enstro- 
phy cascade, i.e. for spatial scales r] u < r < rj n . Separat- 
ing small and large scales one finds energy and enstrophy 
balance equations where terms resulting from the small 
scale fluctuations act like an effective eddy viscosity for 
the large scale components of lo. Analytical expressions 
for the 2nd order vorticity structure function ( r ) an( l 
the 2nd order velocity structure function D 2 {r) can be 
found using the Batchelor interpolation technique |^6],[36| . 
In dimensionless form they read 



(1 + af 2 ) 1 / 3 




D^(oo) 



(14) 
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with the parameter a 



2/3 



■^2 an( i the asymptotic value 
= 14.8. This spectrum also de- 



pends on the energy dissipation e, which when expressed 
in the length and energy scales of the enstrophy cascade 
becomes the dimensionless parameter Re* = e/(et/ 3 u). 
The structure functions are shown in Fig. 1. Besides the 
prominent f 2 behavior that follows already by dimen- 
sional analysis one notes an intermediate scaling with 
f 4 / 3 ; the range over which this scaling is observed de- 
pends on Re* (see below). The corresponding longitudi- 
nal velocity structure function Dn(f) is given with eq. 
©by 
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(1 



r 2 )5/3 



(1 + 0^)2/3 _ l 



bar 2 



2af 2 



(15) 



For the energy ISR no such analytical expression is 
known. We therefore combine a model for the energy 
distribution in fc-space with numerical transformations 
to obtain the longitudinal structure function. Recent ex- 
periments on forced two-dimensional turbulence [p|,p"2"| , a 
number of direct numerical simulations [pOLp7H41fL field 



theoretical investigations |42| as well as cascade models 
p3| support the existence of a Kolmogorov-like scaling 
for the energy spectrum, E(k) ~ fc -5 / 3 for (fc < kf), in 
the energy ISR and E{k) ~ k" 13 with (3 > 3 for (k > kf) 
for the enstrophy ISR. We therefore start with the fol- 
lowing model spectrum for the amplitudes (|uk| 2 ) of the 
velocity field in a Fourier representation in a periodic box 
of size L = 2tt 



(|u k | 2 }~ 



k-P 



exp 



A- 3 

fc -2/3 

2" 



h < k < k f , 

kf <k<k u = JL, /?>2, 



(16) 



Note the different scalings for (|uk| 2 ) and the energy 
spectrum E(k) due to phase space factor, i.e. E(k) ~ 
k^ 13 ^ 1 corresponds to (|uk| 2 ) ^k" 13 . 

The first range approximates finite system size effects 
where we have chosen a slope of 3 in correspondence with 
results of numerical experiments |3^^l| ■ This is followed 
by the inverse energy cascade range with a Kolmogorov- 
like scaling law. At the injection scale kf the enstrophy 
cascade to larger values of k starts, followed by the vis- 
cous cutoff. The energy spectra with (3 = 3 for three dif- 
ferent values of the injection wavenumber kf are shown 
in Fig. |. 

The relation between velocity spectrum scaling and the 
velocity structure function -D2M assuming stationarity, 
homogeneity, and isotropy is given by the volume average 

D 2 (r) = ± J |u(x + r)-u(x)| 2 oV, 

= T7 I lE u " ex P( ik ' x )I ex P(* k ' r ) - 1] I" 
v Jv k 

= 2]T(|u k | 2 )(l-cos(k.r)). (17) 

k 

By averaging over all directions (due to isotropy) in k- 
space the cosine gives rise to the Bessel function Jg(fcr), 



# 2 (r) = 2]T<|u k | 2 )(l-Jo(fcr)). 



(18) 



The model spectrum (|l6| ) is then substituted and the 
summation in (|l^) is evaluated numerically using a fi- 
nite, geometrically scaling set of wave numbers. It should 



be mentioned here that the model does not contain a 
spectral range that would correspond to the intermedi- 
ate r 4 / 3 -scaling of the structure function in the enstrophy 
ISR found in the analytical theory. We will come back to 
this point in the discussion of our results. 



III. RESULTS 



A. Fractal dimension in the enstrophy ISR 



We first calculate the scaling behavior in the enstro- 
phy ISR where the analytical expression (Q) is available. 

Inserting ([li]) in ([l2|) and neglecting the term ^3 for the 

(21 

moment, one notes that 5 g depends on three quantities: 
the parameter Re* , the Prandtl number Pr, and the scale 



f itself. The numerical result for S 



(i) 



?(2) 



1 are shown 



in Fig. |3j for a Prandtl number range varying over ten or- 
ders of magnitude and Re* = 7.6. The grey shaded area 
denotes the range of scales where f 4 / 3 gives the main 
contribution to the structure function. It is only in this 
range that we find 1 < 4 X) < 2. The range is bounded 
by f\ < f < fi where f\ is the crossover scale from the 
viscous subrange (VSR) and ?2 is the crossover scale to 
the r 2 -scaling in the enstrophy ISR, 



n =-L0D^(oa )-',//, 



f 2 = -L(DW(oo))\Re? -D^(cx,)/2)- 3 / 2 , (19) 
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where D^\oo) = 14.8 has to be taken. The larger Re* 
the smaller the range of the f 4 / 3 -scaling. It can be ob- 
served only for Re* within the interval 



7.4? 



< Re* < 



:25. (20) 



The lower bound follows from the positivity of the struc- 
ture function by its definition (cf. second term of (pT[)). 
The upper bound is a result of eq. ( |l^ ) and the constraint 
fx > 1. For Re* approaching 7.4 follows f-i goes to infin- 
ity. The f 4 / 3 - scaling range is then extended over the 
whole enstrophy ISR. We see in Fig. |4| that for increasing 
Re* the intermediate fractal scaling of the graph is more 
and more suppressed and conclude that this behavior of 
Sg is due to the presence of the f 4 / 3 -scaling range. The 
above estimates give f\ ~ 6.0 and f 2 ~ 1400 for Re* = 7.6 
and fi « 1.0 and ?2 ~ 1.7 for Re* = 25.0, respectively. In 
the lower panel the corresponding scaling exponent of the 



„(0) 



4 1} is plotted. The 



scalar structure function Q 

plateau of the structure function for large Prandtl 
number and scales below the smallest scales in the turbu- 
lent fluid irjj]^ < 1) corresponds to the Batchelor regime 
of chaotic scalar advection in a smooth fluid |1| . 

For small values of Pr the diffusion k dominates the 
passive scalar dynamics. The scalar field is smooth, Sg 1 ^ = 



1. The exponent 5 g 
square root of eq. 



j;rows when the second term in the 
becomes dominant. By inserting 



the power law D\\ = ■^(D 2 ul \oo)r) i / 3 ' for the enstrophy 



ISR at 



one gets a crossover for 



r c = i y^Fr- 3 / 5 (^(oo))- 2 / 5 «0.48Pr- 3 / 5 . (21) 

By putting f c = r-i and using dl3) the maximum Prandtl 
number Pr s without fractal 5 g can be estimated as 



The structure function of a passive scalar in the enstro- 
phy ISR shows four different regimes. For very small f 
smoothness gives Sg 1 * 1 = 1. This is followed by the Batch- 
elor regime Sg 1 ^ = 2 for sufficiently large Pr. The f 4 / 3 - 
scaling discovered by Grossmann and Mertens is reflected 
in a decrease of 5g below 2 near r/rj u w 10 lztl . For larger 
f it goes back up to 2. 

So far we neglected the term ^3 = aPrr 2 (see (Q) 
and (0)) in our calculation. Because of its f 2 -scaling it 
dominates the structure function for large f. In this 
term was assumed to be subdominant. Substituting the 
various definitions it can be expressed as a ratio of two 
rates, 



IV0I 



(i/flVwl 2 )) 1 ^) 



TV, 



(25) 



- 1/3 



The rate rg — eg/6 2 ms is a scalar forcing rate. r b 
is the strain rate in the enstrophy cascade and character- 
istic of the passive scalar advection by the vortices. The 
case a > 1 then corresponds to rg > r u , i.e. fast driving 
and slow advection. Then the scalar field fills space and 
~ 2. In the other case, a < 1, the advection domi- 
nates and the structure function of the fluid is reflected in 
that of the scalar. It is this latter case that was discussed 
in [gij for surface waves. The size of a is determined 
by the experimental situation and has to be taken from 
measurements. All quantities that enter ( |25| ) are experi- 
mentally accessible; note that the enstrophy dissipation 
rate is related to velocity gradients via e w = — 8((d x u x ) 3 ) 
& 

Results for different Pr with a = 1 are shown in Fig. 0. 
The main effect of an increasing ^3 is the suppression of 
the crossover scaling and a transition for large r. 



B. Extension to the energy ISR 



Pr s < 2^5(L)f ) (oo))- 4 (i?e* - ^ ) (oo)/2) 



N5/2 



(22) 



With Re* = 7.6 and 25.0 this gives Pr s < 2 ■ 10" 6 and 
10 _1 , respectively. 
For large values of Pr one observes a transition to 
2 even when the velocity field is in the VSR. Again 



the second term of (|12j) dominates because of its large 

Re* 
8 



prefactor Pr. Taking Du = ^§-f 2 for the VSR gives 



With f c = we get those Pri which give S g 
the VSR over at least one decade of scales, 

Pri > 200V6( J D^ ) (cx)))- 4 (i?e*) 5/2 . 



(i) 



(23) 
2 in 

(24) 



For Re* = 7.6 and 25.0 this results in Pr ; > 2.0 and 30.0, 
respectively. 



The extension of D2 (r) to the whole range of scales is 
done with ( pj| ) and the results for Sg 1 ^ are given in Fig. ^ 
for three input model spectra (see Fig. ||) which differ by 
the injection wavenumber kf. The smaller kf the longer 
is the enstrophy ISR extended which results in a domi- 
nant range where 5 g = 2. On the other hand, the larger 
kf the more dominant the inverse energy cascade range, 
indicated as the grey shaded area in Fig. |. The corre- 
sponding longitudinal velocity structure function D\\ (r) 
is superimposed. Note that the model spectrum has to 
be normalized to give D\\ = ^§-r 2 in the VSR. In the 
enstrophy ISR we find Du (r) ~ f 2 and in the energy ISR 

D\\(r) ~ r 2 / 3 , leading to Sg 1 ^ — 2 and Sg — |, respec- 
tively. As mentioned, the model spectrum does not show 
the f 4 / 3 -scaling predicted by [^(|. Therefore, if Re* is in 

the range where a f 4 / 3 -scaling appears the Sg 1 ^ values for 
r ~ r/uj have to be replaced by the ones in Figs, ^, and 
0. For very large values of r we can replace Jo(fcr) by 



5 



^ cos(fcr - f) resulting 

in D 2 (r) ps 2 (|u k | 2 ) in (|l8|). The constant asymp- 
totic behavior of the structure function corresponds with 



its asymptotic form J (kr 

|2\ 



(Cf. ©). 



The model spectrum contains a free parameter (3 which 
has no agreed upon value. Numerical simulations |38] Q 

suggest a range /3 € [2,4]. For /? = 2 we get £g slightly 
below 2 in the enstrophy ISR which changes clearly to 
4^ = 2 for /3 > 2 (cf. Fig. 0). As expected, the value of 
5g in the energy ISR is insensitive to a /3-variation. 

Again we have to discuss the additional influence of 
the 73 term in ([?]). Will inverse cascade effects be sup- 
pressed in the large Pr number case because of the dom- 
inance of the f 2 -scaling at large separations? In order 
to determine the scale f c where I3 > 1%, we use the ex- 
perimental value for the Kolmogorov constant Ck || and 
assume a completely extended inverse cascade with no in- 
termittency corrections. Then D 2 (r) = 4Ck e 2 / 3 / °°(1 — 
J (fcr))fc- 5 / 3 dfc and D 2 = 6 2 r 2/3 . With C K between 5.5 
and 7, we find b 2 between 31.5 and 40 for the energy ISR 
and thus finally 



> 



8a 2 ' 



(26) 



where I lies between 14 and 17. The scale f a is shifted 
towards larger values for decreasing a. A factor a~ 1 can 
suppress the scaling behavior in the energy ISR which 
was found above completely. This fact is illustrated in 
Fig. ||. Clearly the asymptotic state for f to infinity leads 
here to Sg 1 ^ approaching 2. 



IV. DISCUSSION 

Our main findings for a passive scalar in a 2-d turbu- 
lent flow field can be summarized as follows: (1) There is 
a critical scale set by equation below which the spec- 
trum is smooth, Sg = 1, because of diffusion dominance. 
(2) Between this scale and the injection scale ri n the scal- 
ing exponent Jg 1 ' =2 in most cases. (3) An exception is 
found for Re* in the interval set by (EG) , where a scaling 

exponent 5g < 2 is found. The limits of this interval are 
given by ( |19|) and the deviation from 2 is controlled by 
the parameter a, eq. (p5|). (4) Beyond the injection length 
and up to a length set by equation (^), the scalar field 
scales with the exponent Sg^ = 5/3 as expected for the 
energy inertial subrange. (5) Above the length scale set 
by (|2^), the exponent again increases to 2. What is most 
surprising is that the scaling derived within geometric 
measure theory depends not only on the scaling of the 
velocity field but also on two additional dimensionless 
numbers, the Reynolds number Re* which causes the in- 
termediate scaling in the enstrophy viscous subrange and 
on a which suppresses the velocity field induced scaling 
at large separations for rapid driving. 



At this point input from experiments on two- 
dimensional turbulence is necessary to check and expand 
the theoretical results. Cardoso et al. || measured disper- 
sion in a quasi-two-dimensional turbulent flow and com- 
pared with results for the energy inertial subrange. They 
observed a velocity structure function with scaling r° and 
a fractal dimension 5 g X ^ between 1.3 and 1.5 with an aver- 
age of about 1.4. Substituting a velocity scaling function 
Dm = Cf° in our main equation (fi"^) gives 
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(27) 



If the quadratic term can be neglected, i.e. if f is small 
enough, the inequality reads Sg 1 ^ < 3/2. The experimen- 
tal results are indeed below but close to this limit, so that 
the assumption that the distances are small is probably 
reasonable. For larger separation there is a crossover to 
^g 1 "* < 2, and it would be interesting to see whether the 
experimental data follow this behavior. For the energy in- 
ertial subrange [and not too large separations, see (p6|)], 
the inequality would be 5 g < 5/3, higher than the one 
for the experimentally observed spectrum. 

Further experiments or numerical studies to check the 
results from geometric measure theory, especially the 
ones for the enstrophy cascade and for the dependence 
on a, are clearly needed. Perhaps it is possible to com- 
bine the experiments on passive scalar mixing with 
the set-up for extended, stationary inverse and direct cas- 
cades p|,|3"H in order to measure the scaling behavior men- 
tioned in (2). In order to check the predictions for the 
enstrophy cascade in (1) the spatial resolution has to be 
enlarged. Otherwise e.g. the existence of the intermediate 
f 4 / 3 -scaling of D 2 if) cannot be detected. We remind the 
reader that this range is only well established for values 
of Re* close to its lower threshold (see Fig. Q). Its local- 
ization with respect to f prevents it from being seen in 
the Fourier spectrum, as already discussed by Grossmann 
and Mertens p6[ . 

Another open question which calls for more input from 
numerical simulations and experiments is that of the 
scalar flatness in 2-d. For a non-Gaussian scalar statistics 
we would expect a scale-dependent flatness Fg(r) caus- 
ing a further scale dependence of the third term in ( |l^ ) 
and thus leading to a modification of the present model. 

The problem studied here has also interesting links to 
magnetohydrodynamics. First steps towards using geo- 
metric measure theory in this context were undertaken 
by Grauer and Marliani (44| . In two dimensions there is a 
direct relation between magnetic field advection and the 
scalar dynamics studied here since the vector potential 
for the magnetic field has only a z-component. Conse- 
quences of this relation are under investigation. 



6 



[1] S. Edouard, B. Legras, F. Levevre, and R. Eymard, Na- 
ture 384, 444 (1996). 

[2] D. K. Lilly, J. Atmos. Sci. 46, 2026 (1989). 

[3] M. Lesieur, Turbulence in Fluids, (Martinus Nijhofl Pub- 
lishers, Dodrecht) 1987. 

[4] J. Sommeria, J. Fluid Mech. 170, 139 (1986). 

[5] P.Tabeling, S. Burkhart, O. Cardoso, and H. Willaime, 
Phys. Rev. Lett. 67, 3772 (1991). 

[6] J. Paret and P.Tabeling, Phys. Rev. Lett. 79, 4162 
(1997). 

[7] B. S. Williams, D. Marteau, and J. P. Gollub, Phys. Flu- 
ids 9, 2061 (1997). 
[8] O. Cardoso, B. Gluckmann, O. Parcollet, and P.Tabeling, 

Phys. Fluids 8, 209 (1996). 
[9] M. Gharib and P. Derango, Physica D 37, 406 (1989). 

[10] B. K. Martin, X. L. Wu, W. I. Goldburg, and M. A. Rut- 
gers, Phys. Rev. Lett. 80, 3964 (1998). 

[11] M. Rivera, P. Vorobieff, and R. E. Ecke, Phys. Rev. Lett. 
81, 1417 (1998). 

[12] M. A. Rutgers, Phys. Rev. Lett. 81, 2244 (1998). 

[13] R. H. Kraichnan, Phys. Fluids 10, 1417 (1967). 

[14] G. K. Batchelor, Phys. Fluids Supplement 2, 233 (1969). 

[15] G. K. Batchelor, J. Fluid Mech. 5, 113 (1959). 

[16] A. M. Obukhov, Izv. Akad. Nauk SSSR, Ser. Geog. Gc- 
ofiz. 13, 58 (1949). 

[17] S. Corrsin, J. Appl. Phys. 22, 469 (1951). 

[18] R. Benzi, S. Patarnello, and P. Santangelo, Europhys. 
Lett. 3, 811 (1987). 

[19] A. Babiano, C. Basdevant, B. Legras, and R. Sadourny, 
J. Fluid Mech. 183, 379 (1987). 

[20] C. Basdevant and T. Philipovitch, Physica D 37, 17 
(1994). 

[21] P. Constantin, I. Procaccia, and K. R. Sreenivasan, Phys. 

Rev. Lett. 67, 1739 (1991). 
[22] P. Constantin and I. Procaccia, Phys. Rev. E 47, 3307 

(1993) . 

[23] P. Constantin and I. Procaccia, Nonlinearity 7, 1045 

(1994) . 

[24] I. Procaccia and P. Constantin, Europhys. Lett. 22, 689 

(1993) . 

[25] S. Grossmann and D. Lohse, Europhys. Lett. 27, 347 

(1994) . 

[26] S. Grossmann and P. Mertens, Z. Phys. B 88, 105 (1992). 

[27] K. J. Falconer, The Geometry of Fractal Sets, (Cam- 
bridge University Press, Cambridge) 1985. 

[28] H. Federer, Geometric Measure Theory, (Springer, 
Berlin) 1969. 

[29] F. Morgan, Geometric Measure Theory, a Beginners 

Guide, (Academic Press, Boston) 1988. 
[30] L. M. Smith and V. Yakhot, Phys. Rev. Lett. 71, 352 

(1993) . 

[31] J. Paret, M.-C. Jullien, and P.Tabeling, Vorticity Statis- 
tics in the two-dimensional enstrophy cascade, submitted 
to Phys. Rev. Lett. (1999). 

[32] B. I. Shraiman and E. D. Siggia, Phys. Rev. E 49, 2912 

(1994) . 

[33] R. H. Kraichnan, Phys. Rev. Lett. 72, 1016 (1994). 
[34] A. S. Monin and A. M. Yaglom, Statistical Fluid Mechan- 
ics, (MIT Press, Cambridge, Massachusetts) 1975. 
[35] L. D. Landau and E. M. Lifschitz, Course of Theoretical 



Physics Vol. 6, (Pergamon Press, Oxford) 1987. 
[36] G. K. Batchelor, Proc. Camb. Phil. Soc. 47, 359 (1951). 
[37] U. Frisch and P. L. Sulem, Phys. Fluids 27, 1921 (1984). 
[38] R. Benzi, C. Paladin, S. Patarnello, P. Santangelo, and 

A. Vulpiani, J. Phys. A 19, 3771 (1986). 
[39] V. Borue, Phys. Rev. Lett. 71, 3967 (1993). 
[40] N. K.-R. Kevlahan and M. Farge, J. Fluid Mech. 346, 49 

(1997). 

[41] A. Babiano, B. Dubrulle, and P. Frick, Phys. Rev. E 55, 
2693 (1997). 

[42] G. Falkovich and V. Lebedev, Phys. Rev. E 49, R1800 
(1994). 

[43] J. Schumacher, Diploma thesis, Philipps University Mar- 
burg, 1994 (unpublished). 
[44] R. Grauer and C. Marliani, Phys. Plasmas 2, 41 (1995). 





8 










/ 




6 


^VSR 




ISR 


/' ' 




CV 3 
Uj 
W 3 








// .■■ 




P- 


4 






/' s 






m 
o 


2 



-2 


/ V 
'/ 


/// 




Re*- 
--- Re*- 
Re*- 
Re"- 


7.4 
7.6 
1 1.0 
25.0 



-10 12 3 4 
log 10 (r/7jj 

FIG. 1. Velocity structure function D2(r) in the enstrophy 
inertial subrange for four different values of Re*. 
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FIG. 2. Model spectrum (|u k | ) for three different values of 
k f indicated by the arrows (k f « 2 • 10~ 5 , 5 ■ 10~ 4 , 10~ 2 ). The 
wavenumbers are given in units of k^ — 77J 1 . The exponent /3 
was set to three. 
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FIG. 3. Fractal dimension Sg 1 ^ for three Prandtl num- 
bers and the corresponding velocity structure function D2(f) 
(thick line) for Re* = 7.6. The grey shaded area denotes the 
range of scales where the f 4//3 term dominates for the param- 
eter set. A fractal Sg 1 ^ can be observed in this range of scales. 
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FIG. 4. Fractal dimension & g X) and scaling exponent C}^ 
as a function of Pr = 10 A and of Re*. The solid line is 
Re* = 7.4 (the lower bound), the dotted line is Re* — 7.6, 
and the dashed line is Re* — 13.0 in both panels. Note that 
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FIG. 5. Fractal dimension Sg 1 ^ as a function of Pr = 10 A 
for Re* = 7.6, a = 1. The solid line plots show the results 
when only the advection term Ii is taken. The dotted lines 
show the additional influence of the forcing term I3. 
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FIG. 6. Longitudinal velocity structure function D\\ (f) 
(thick line) and fractal dimension Sg 1 ^ over f for Re* — 7.6 
and Pr = 10 4 , 10 2 , 10°, 1(T 2 , and 1(T 4 decreasing from left 
to right. The grey shaded range of scales denotes the in- 
verse energy cascade range of (|u k | ). (a): kf w 2 • 10 , (b): 
kf ?s 5 ■ 10~ 4 , (c): kf « 10~ 2 . The exponent (3 = 3 was taken. 
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FIG. 7. Fractal dimension 5^ for Pr = 10 4 and Re* = 7.6 FIG. 8. Fractal dimension for Re* = 7.6 and for three 

for different values of the scaling exponent (3 taken in the different values of the parameter a — ege^ 1 ^ 3 /Or.m.s.- The 
enstrophy ISR for the model spectrum (|uk| 2 ) (cf. (Q). exponent (3 = 3 was taken. 
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